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Abstract We report and analyze the observational evidence of global kink oscillations in a 
solar filament as observed in Ha by National Solar Observatory (VS’Oj/Global Oscillation 
Network Group (GONG) instrument. An Ml. 1-class flare in active region 11692 on 2013 
March 15 induced a global kink mode in the filament lying in the south-west of AR11692. 
We find periods of about 61-67 minutes and damping times of 92-117 minutes at three 
vertical slice positions chosen in and around the filament apex. We find that the waves are 
damped. From the observed global kink mode period and damping time scale using the 
theory of resonant absorption we perform prominence seismology. We estimate a lower 
cut-off value for the inhomogeneity length-scale to be around 0.34—0.44 times the radius 
of the filament cross-section. 
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1 INTRODUCTION 


Solar prominences manifest various kinds of magnetohydrodynamic (MHD) oscillatory motions 
(Arregui et al. 201 2) . These oscillations are broadly classified either as large amplitude oscillations 
where prominence as a whole oscillates with a velocity amplitude of a few tens of km s _1 or as small 
amplitude oscillations localized in a part of a prominence and with a velocity amplitude of about 2-3 km 


s 1 . Many cases of prominence thread oscillations with small amplitudes have also been reported (Yi 


et al. 


1991 Yi & Engvold | 1991). Horizontal flows with simultaneous transverse oscillations were re¬ 
ported by Lin et al. ( 2002) , |Lin et al. ( 2005) and Okamoto et al. ( 2007) . Okamoto et al. ( 2007) reported 
the propagating Alfven waves in the fila ment th reads and calculated the wave velocity to be about 1050 


km s . However, later Terradas et al. (2008) correctly interpreted them as standing kink waves and 


Lin 


calculated the lower bound of Alfven velocity in these threads with a lower value of 120 km s 1 
et al. (2009) also reported the propagating kink waves in the filament threads. Ning et al. (2009) found 


vertical and horizontal oscillations in the prominence threads with simultaneous drifting in the plane of 
sky. 
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There have been numerous reports on small amplitude oscillations as summerized in a review article 
by Oliver & Ballester (2002) and large amplitude oscillations (Tripathi et al., 2009). Ramsey & Smith 
[ ( 1966) reported the observational evidence of large amplitude oscillations induced by disturbances 
coming from a nearby flare. Recent observations of large amplitude oscillations were reported in Ha 


(Eto et al. 2002 Jing et al. 2003 Okamoto et al. 2004) Jing et al. 2006 Liu et al. 201 3[>, in E xtreme 


Ultraviolet (Isobe & Tripathi 2006)) and in He 10830 infrared emissions (Gilbert et al. |2008|. Most 

recent report was by Hershaw et al. ( 201 1) in which they investigated two large amplitude damped 
transverse oscillations in an EUV prominence on the solar bmb. The large amplitude oscillations can be 
either longitudinal (Jing et al. 2003) 2006) Vrsnak et al. 2007 Li & Zhang |2012| Luna & Karpen 


2012 Luna et al. 2012 Zhang et al. 2012 Luna et al. 2014) or transverse (Isobe & Tripathi 2006 
Okamoto et al. 2004 Schmieder et al. 2013). Detection of magnetohydrodynamic (MHD) waves and 


oscillations provide important input in diagnosing the local plasma conditions of the solar corona by 
applying the principle of MHD seismology (Nakariakov & Verwichte 2005 Andries et al. 2005, 


2009) |Ruderman & Erdelyi | 2009). The fundamental mode of the global kink wave is characterized 


by a displacement of the magnetic flux tube where all its parts are in phase (Nakariakov et al. 1999). 
In the MHD regime, the global transverse oscillations of magnetic flux tubes are interpreted as nearly 
incompressible fast kink modes (Edwin & Roberts , 1983 Roberts et al. 1984) Erdelyi & Taroyan") 
|2008) [Van Doorsselaere et «/.||2008 Goossens et al. |2009). 


Prominence oscillations are often found to be damped (Ruderman & Roberts 2002| Erdelyi & 


Taroyan 1 12008 Andries et al. 2009) |Arregui & Ballester | 201 1) . Damping is observed in both large 
amplitude oscillations where the prominence oscillates as a whole (Hershaw et al. 2011) and at smaller 


2002 


scales, where different parts of the prominence show different damping timescales (Terradas et al. 

|Lin et al\ |2009) >. 

Besides the oscillation period, the damping timescale also provides fundamental information about 
the physical conditions around the prominence. The MHD seismology of the localized plasma environ¬ 
ment of the magnetic flux tubes based on both the oscillation period and damping timescale of transverse 
waves is performed in a consistent manner by Goossens et a/.|( 2006); |Arregui et al. (2007); Goossens 


[ (2008); Goossens et al. (2008); Soler et al. (2013). These fast kink radial modes get attenuated likely 


due to the resonant absorption and are converted into dominantly torsional (azimuthal) highly incom- 


pressible motion (P 

alderman & Roberts 2002 Goossens et al. 2002, 2006, 2008 2013|i Goossens 

2008; Arregui et al. 

2008). However, other mechanisms (e.g., dissipation through wave leakage, curva- 


ture effects, phase-mixing, etc.) can also be at work in dissipating such waves in the solar coronal tubes 
(e.g, Ruderman & Roberts 2002) Ofman & Wang |2002[ Srivastava et al. 2013) . 


In this paper, we present evidence of large amplitude, long period damped transverse oscillations 
in a Solar filament as observed by GONG instrument of National Solar Observatory (NSO) in Ha. We 
interpret them as global kink waves as they displace the prominence tube as a whole in the transverse 
mode. Such an analysis has been carried out by Hershaw et al. ( |201l) , but they have reported the 
kink oscillations in EUV and there were differences in the periods of the two legs of the prominence. 
With that, the two legs of the prominence oscillated in phase initially and then gradually moved out of 
phase. Okamoto et al. (j2004[) and Eto et al. (| 2002) observed winking filaments in Ha with intermediate 
periods (20-40 min). However, there was no signature of damping. Isobe & Tripathi (2006 ) and |Chen] 


et al. (2008) also reported the undamped transverse oscillations, which happened during or immediately 


prior to eruption. Gilbert et q/.| ( |2008 1 reported long period, large amplitude vertical oscillations but the 
damping time observed by them is six times the period of the oscillation. Prominence seismology in 


large amplitude oscillations has been carried out by Pinter et al. (2008) and Vrsnak et al. (2007 1 , and 
in prominence threads by Lin et al. (2009 ). The paper is organized as follows: In Sect. 2, we present 
briefly the observations. We describe the detection of kink oscillations in Sect. 3. In Sect. 4, we perform 
the MHD seismology and estimate the localized physical conditions within the filament. In Sect. 5, the 
discussions and conclusions are outlined. 













































































































































































































































MHD Seismology of a loop-like filament tube by observed kink waves 


3 


2 OBSERVATIONS 


An Ml.1-class flare is observed by GOES satellite in the active region AR11692 on 2013 March 15, 
which is associated with a halo coronal mass ejection (CME). Soon after the onset of the flare, oscilla¬ 
tions are observed in the filament lying in south-west to the active region as seen in Ha from NSO/GONG 
(Fig. □ bottom panel, the movie is available online. 

NSO/GONG provides full-disk observations of the Sun at 6563 A. It has a maximum pixel resolution of 
~ 1.07" and a cadence of 1 minute. The flare or the associated CME generated global disturbances that 
in turn triggered the observed transverse oscillations in the body of the filament. Observations started 
at 06:00:54 UT and ended at 10:22:54 UT. At least four transversal cycles of the global transverse os¬ 
cillations are observed with their significant damping. The top panel of Fig.[T]shows the Helioseismic 
and Magnetic Imager (HMI) magnetogram contours overlayed on the image observed by Atmospheric 
Imaging Assembly (AIA) on board SDO at 304 A. AIA observes the Sun in seven Extreme Ultraviolet 
(EUV) bandpasses and has a pixel size of 0.6"(Lemen et al. 2012). The 304 A images are chosen 
because in this passband the footpoints of the filament can be identified. Footpoints are identified by 
examining the strands of the filament in the image together with contours of opposite polarities because 
a filament always lies on the polarity inversion line. Footpoints are marked as ‘X’ (negative) and ‘Y’ 
(positive) on the image. Green and black contours represent constant magnetic field strength of 20 G 
and -20 G, respectively. We chose five points over the filament between the foot points ‘X’ & ‘ Y’. We 
then interpolate the curve using cubic spline between the footpoints. The interpolated curve has to pass 
through the chosen five points. The length of the interpolated curve should be approximately equal to 
the length of the filament. We repeat this process several times by choosing different five points every 
time between ‘X’ and ‘Y\ Therefore, the length of the interpolated curve will vary, thus its mean value 
and the standard deviation are estimated. The length of the filament, L , is measured to be ~ 235±8 Mm. 
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Fig. 1 The left panel shows the Helioseismic and Magnetic Imager (HMI) magnetogram 
contours overlayed on the image as observed by Atmospheric Imaging Assembly (AIA) on 
board SDO at 304 A channel. The right panel shows the NSO/GONG Ha image of the 
filament south-west to AR 11692. Slices used for creating time-distance map are marked as 
white, yellow and black vertical lines respectively. 


3 DETECTION OF GLOBAL KINK OSCILLATIONS IN A FILAMENT TUBE 

The evolution of the oscillations is obtained by placing three artificial slices close to the apex of the 
filament parallel to the observed oscillations and perpendicular to the filament (see the right panel of 
Fig. 0 using the NSO/GONG Ha image sequence. For each slice a two dimensional time-distance 
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diagram is created (see the left panel of Fig. [2]», where the x-axis represents the time in minutes and the 
y -axis represents the distance along the slice in Mm. Thick black region represents the data gap from 
07:07:54 UT to 07:18:54 UT. Along each column of the time-distance diagram a Gaussian curve is fitted 
and the mean values with one-sigma error bars are estimated. The entire time-distance diagram is then 
fitted with a damped harmonic curve represented by the formula (Aschwanden et al. 1999|) 


y(t) = c + aosin(ojt + <f>)e t ^ Td , 


( 1 ) 


where c is a constant, ao is the amplitude, r,/ is the damping time, ui is the angular frequency, and d is 
the initial phase. 

The least square fitting is performed using the function MPFIT.pro ( |Markwardt 2009]) in the 
Interactive Data Language (IDL). The best-fit curve is shown in Fig. [2] (right panel). The period of 
the oscillation P = —, the damping time r,/, and the initial phase <p for the first (second and third) 
slices are found to be 67±3 min (63±2 min and 61 ±2 min), 98 ±47 min (117 ±45 min and 92±28 min) 
and 2.7±0.4 (2.9±0.3 and 2.7±0.3), respectively. From Fig. [2] (right panel) we note that the filament 
segments at the three slices are oscillating in phase. The amplitude of the oscillation, ao, is 8±4 Mm 
(9±3 Mm and 13±3 Mm) and the velocity amplitude is about 12±6 km s” 1 (15±6 km s” 1 and 22±7 
km s” 1 ) for the first (second and third) slice, respectively. Oscillations with decaying amplitudes are 
also noticed near the footpoints of the filament. They are in phase with the oscillations seen near the 
apex of the filament. The in-phase displacement of the filament is interpreted as the global standing kink 
mode. Since all parts of the filament oscillate in phase and the footpoints of the filament are fixed, thus 
we infer it to be the fundamental standing mode. Moreover, it should be noted that we can not distin¬ 
guish on the observational basis whether the oscillations are vertical or horizontal, or mixed. However, 
this does not affect the principle of the MHD seismology of the filament that we intend to implement 
in the present paper. Assuming the filament to be a flux tube embedded in uniform magnetic field, the 
phase speed for the fundamental standing mode can be calculated using the relation. 


Vph — 


2 L 
~P' 


( 2 ) 


As a result, v p h for the first (second and third) slice is computed to be 117±9 km s -1 (124±8 km s -1 
and 128±9 km s -1 ). The sound speed, c s , with the chromospheric temperature is 15 km s -1 . With 
Vph > c s , we infer it to be the fast kink oscillation. 

Prominence models in general invoke the fine structure of isolated/groups of prominence threads on 
which the MHD wave mode(s) distribute ( |Arregui & BallesteF 20111. In the present case, almost bipolar 
prominence having the aspect ratio rather small, like the EUV loops, is subjected to the transverse 
displacement of its axis as a whole. Therefore, we treat it as a loop-like prominence tube supporting the 
kink oscillations. 


4 CORONAL SEISMOLOGY OF THE OSCILLATING FILAMENT TUBE 


We assume the filament to be a flux tube with its radius much smaller than its length, L (i.e., thin-tube 
approximation). With the thin tube approximation or long wavelength limit (i.e., k z a « 1, where 
k z = 2 ? and a is the cross-section of the filament, for the fundamental mode, A=2L), the kink wave 
speed ( Vk ) is approximately equal to the phase speed (Vph)- 

Analytical relation between the period of oscillation P, the internal Alfven travel time t a ^ and the 
density contrast Q is given by 


_ TA,i __ 1 , C 
y p V2 ( C±1 


) 1/2 , 


(3) 


where t a ,i = L/vaa, the internal Alfven speed VA,i = F is the permeability of the medium, and 


Pi is the density inside the filament (Goossens et al. 2008). We take p = Air X 10 7 H??r 1 in SI units. 
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Fig. 2 Left, time-distance maps (top to bottom) corresponding to the three slice locations as 
marked in Fig. 1 by white, yellow and black lines respectively. Right. Exponentially damped 
sinusoidal curves 
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£ in Eq. (3) is the internal to external density ratio —. From Eq. (3), the Alfven speed inside the filament 
tube is given by 


«a,. = 


(4) 


The damping time of the fast kink waves may give an insight of ambient plasma and structure of 
the filament. [Arregui et q/.| ( |2008) proposed resonant absorption as a damping mechanism in the context 
of prominence fine structure oscillations. With the thin boundary (TB) approximation, the analytical 
expression of the damping time in the asymptotic limit in the Cartesian coordinates was derived by 
Hollweg & Yang (1988]) and later by Goossens et al. ( 1 992) and Ruderman & Roberts (2002]). The 
asymptotic expression is (see|Goossens et al. 2006[ Arregui et al\ |2008|) 


Td 


2 R Pi 


P 7T l Pi — p e ’ 


(5) 


where R is the mean radius of the filament and l is the thickness of its non-uniform layer. Sinusoidal 


variation of the density is assumed across its non-uniform layer leading to the factor of — (Ruderman & 
2002 1 instead of 4 2 for linear variations (Goossens et al. 1992 1 . Using Eq. (5), we can express 


Roberts 


the inhomogeneity length scale (z = 1/2R) as (Goossens et al. 2008) 


z = 


1 C + 1 
CC- 1’ 


( 6 ) 


where C is . 


However, the first paper that consistently used information on period and damping time is by Arregui 


et al. (2007). It used the full numerical results of the eigenvalue computations and its scheme is fully 
numerical. From Eqs. (4) and (6), we note that there are three unknown parameters va . i . C> and z. There 
are infinite choices of va, i, (, and z available to satisfy Eqs. (4) and (6). Therefore, we have to estimate 
one unknown and must express the other two unknowns in terms of it. We consider to fix the density 
contrast (Q and express the Alfven speed (va.i) and inhomogeneity length scale ( z ) in terms of £. 

Maximum value of ( is infinity. Since z is the decreasing function of £, z attains the lowest possible 
value for the largest possible value of (. Thus z rnlrl = ^ for C, —>■ oo. The possible maximum value 
of z is 1, therefore, the minimum value of £ will be when z= 1, which can be computed using Eq. (6) 

as Cmin = and using Eq. (4) as VA.imax = t ^ le next ste P’ we com P ute the magnetic 

field strength and inhomogeneity length scale separately corresponding to slices I, II and III for a fixed 
value of (. 


4.1 Density estimates using automated DEM analysis 

It is evident that the Alfven speed in a medium depends on magnetic field strength and density of the 
medium. Using Eq. (4), Alfven speed can be estimated for a given value of (. Therefore, magnetic field 
inside the filament can be estimated if density of the filament is known. To estimate density values we 
follow automated temperature and emission measure analysis technique as developed by Aschwanden 
et al. (2013) to derive the average density and temperature inside the filament tube. Using this technique 


8.60±0.32 


cm 


-3 


the electron density (n e ) and temperature (T e ) inside the filament is estimated to be 10 
and io 5 - 93±0 - 25 K (see Fig. 3). Note that the density is underestimated as compared to the typical density 
range of prominences which ranges from 10 9 to 10 11 cm -3 for quiescent and active filaments (Labrosse 


et al. 2010). We feel that because of the low density contrast with respect to the background as visible 


in various SDO/AIA channels, the DEM method does not work very well. Furthermore, the observed 
filament is made-up of mostly cooler material, whereas various AIA channels are sensitive to hotter 
plasmas, thus the contribution from the prominence material would be minimal. The bright points seen 
in AIA 171 A channel (see Fig. 3) correspond to the prominence-corona transition region (PCTR). 
Thus the density values calculated by the DEM technique is the density of the PCTR associated with 
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Table 1 Seismic parameters for exponential damping 


c 

Slice I 

VA,i z 

Slice II 

VA,i Z 

Slice III 

VA,i Z 

1.414 

- 

- 

115 

1 

— 

1.534 

- 

- 

113 

0.814 

116.69 1 

1.557 

106 

1 

112.67 

0.787 

116.37 0.968 

2 

101.25 

0.654 

107.68 

0.515 

111.21 0.633 

3 

95.46 

0.436 

101.52 

0.343 

104.85 0.422 

4 

92.42 

0.363 

98.30 

0.286 

101.52 0.352 

5 

90.56 

0.327 

96.31 

0.257 

99.47 0.316 

6 

89.29 

0.305 

94.97 

0.240 

98.07 0.295 

7 

88.38 

0.290 

94 

0.229 

97.07 0.281 

8 

87.69 

0.280 

93.25 

0.220 

96.31 0.271 

9 

87.14 

0.272 

92.68 

0.214 

95.72 0.264 

10 

86.71 

0.266 

92.21 

0.210 

95.24 0.258 

50 

83.49 

0.227 

88.79 

0.179 

91.71 0.220 

100 

83.08 

0.222 

88.36 

0.175 

91.26 0.215 

1000 

82.71 

0.218 

87.97 

0.172 

90.85 0.211 

OO 

82.67 

0.218 

87.92 

0.171 

90.80 0.211 


this prominence, rather than the bulk of the prominence material. Thus the DEM forward fitting will 
not be accurate enough for the prominence. By DEM analysis we wish to estimate the PCTR magnetic 
field. If we use these density estimates (from DEM) and assuming the number density ratio H:He to 
be 10:1 inside the filament, the magnetic field strength ( B ) for density contrast ( =100 are estimated 
as 0.86±0.21 G (0.91±0.20 G and 0.94±0.20 G) for first (second and third) slice respectively. These 
magnetic field strengths are much lower than as reported in [Mackay et al. ( 2010[ > for the quiescent 
filaments (3-15 G). It implies that the magnetic field of overlying PCTR is lesser compared to typical 
magnetic field of filament. 


4.2 Magnetic field and inhomogeneity length estimation corresponding to Slice I 

Since automated DEM analysis underestimate the density values thus we take typical density values in 
the filament channel which ranges from 10 9 to 10 11 cm” 3 . Using the observed period P= 67±2 min 
as well as the damping time r ( i=98±47 min, we get ^f=1.46±0.76. For simplicity, henceforth in the 
paper we consider the absolute values of the parameters without error bars. For example, we derive the 
extreme values of the seismic quantities as C=4.59, Cmin=l-557, z m i n = 0.218 and VA,i,max=[06 km 
s” 1 . Table [ljshows the value of v a .i and :: for different values of Q. There are infinite choices as shown 
in Table 1, but we take typical value of £ for the filament to be 100. For £=100 the Alfven velocity and 
z = ^ are 83.08 km s ” 1 and 0.222, respectively, thus l/R= 0.444. For typical density values in the 
filament channel (10 9 to 10 11 cm” 3 ), we deduce the magnetic field value to be ~ 1 .36 G, ~4.28 G, and 
~ 13.6 G, respectively. 


4.3 Magnetic field and inhomogeneity length estimation corresponding to Slice II 

Using the observed period P= 63±2 min and damping time-scale r,/= 117±45 min, we get 
p-=1.86±0.77. The seismic quantities C, C, m in, z m i n and VA,i,max are calculated as 5.83, 1.414, 0.171, 
and 115 km s” 1 , respectively. For £=100, the Alfven velocity and z = 57 ? are 88.36 km s ” 1 and 0.175, 
respectively, thus l/R=Q350. Assuming similar densities values as those for slice I, i.e., n e = 10 9 , 10 10 
and 10 11 cm” 3 , the magnetic field strength B for the given density contrast £ =100 is estimated to be 
~1.44 G, ~4.56 G, and ~ 14.46 G, respectively. 
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Fig. 3 Automated DEM analyses measurements : (i) the selected segment of the tube (top- 
right) which is also shown in various channels of the AIA (bottom-left), (ii) The forward 
DEM fitting over the emission derived along the chosen path (top-left), and derived density 
as well as temperature (blocks in bottom-right multipanel plot). 
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4.4 Magnetic field and inhomogeneity length estimation corresponding to Slice III 

Using the observed period P=61±2 min and damping time t,/= 92±28 min, we get T j[ = 1.51 ±0.51. The 
seismic quantities C, C, m in , z m i n , and VA,i,max are calculated as 4.74, 1.534, 0.211, and 116.69 km 
s -1 , respectively. For £=100 the Alfven velocity and z = ^ are 91.26 km s' 1 and 0.215, respectively, 
thus l/R= 0.430. Assuming similar densities as those for slice I, i.e., n e = 10 9 , 10 10 and 10 11 cm -3 , the 
magnetic field strength B for the given density contrast £ =100 is estimated to be ~1.48 G, ~ 4.70 G, 
and ~ 14.85 G, respectively. 

4.5 Magnetic field and inhomogeneity length estimation for infinite density contrast 

Seismology is performed by using both P and td and the two coupled equations (Eqs. 4-6). We keep the 
one unknown, i.e, density contrast fixed, and express the Alfven speed and inhomogeneity length scale 
in terms of it. In the case of prominence plasmas, if we consider the density contrast is large and tending 
towards infinity, we get two uncoupled equations in which one relates the period and Alfven speed, 
while other the damping ratio with the inhomogeneity length scale (l/R). Therefore, the inversion what 
we performed earlier is independent of the density contrast. This is clearly evident from Table 1 which 
shows that various parameters are having almost similar values for higher values of £ (=1000, oc). For 
£=oc for slice I, II, III, the Alfven velocities and ^ ^ are (82.67, 87.92, 90.8) km s _1 and (0.218, 

0.171, 0.211), respectively, thus l/R=( 0.436, 0.342, 0.422). The magnetic field strength estimated for 
£=oo is similar to what is estimated by assuming £=100 to first place of decimal. 


5 SUMMARY AND DISCUSSION 


A dynamical study of the cool filament observed on-disk in Ha has been performed. We report global 
kink oscillations in the filament system as observed in Ha as the whole filament tube oscillates in the 
transverse mode. The oscillation is most likely excited by the large-scale disturbances originating from 
a nearby M 1.1 -class flare. However, we do not aim to understand about the nature of the driver and its 
properties, it is out of the scope of our present article. Time series analysis revealed a period between 
61 to 67 min and a damping time between 92 to 117 min associated with the kink oscillations around 
the apex of the filament system. The phase speed ( v p h ) is found to be between ~ 117 km s _1 to 128 
km s _1 . These are long-period damped fast kink oscillations, as the phase speed is high compared to 
the chromospheric sound speed. We also find that oscillations are damped. Lin et al. ( 2009) has carried 
out prominence seismology in prominence threads and estimated the magnetic field of about 0.1-20 G. 
However, We have carried out consistent seismology using the values of the periods and the damping 
times at the same time. It is worthwhile to emphasize the use of seismology using the values of the 
periods and the damping times as it gives an order of estimate of upper limit of Alfven speed inside the 
prominence (Goossens et «/.~[ [2008) and hence on the estimation of the magnetic field strength. Pinter et 
al. ( 2008| carried out the prominence seismology and calculated the axial component of the magnetic 
field in a polar crown prominence to be 1-5 G. They used twisted flux rope model to calculate the 
magnetic field values. To the best of our knowledge, there have not been any reports yet on consistent 
seismology in large amplitude prominence oscillations, which we implement in the present work to 
diagnose the local plasma conditions of the observed filament tube. The main result from this study 
is the estimation of the inhomogeneity length scale. Assuming a typical density ratio £ of filament to 
be 100, the inhomogeneity length scale, l/R is 0.444 (0.350 and 0.430) for the first (second and third) 
slice, respectively, which suggests that the inhomogeneous layer is quite thick near the apex of the 
filament. The density ratio estimate is not very accurate, since there are infinite choices, so from the 
available information we can say that l/R> 0.436 (0.342 and 0.422) for the first (second and third) 
slice, respectively. We find the magnetic field strength, using the typical density values, in the filament 
to be in the range of ~1 G to 15 G, which is close to what has been reported by Lin et al. ( 2009) > 
for the prominence threads. To best of our knowledge, this is the first time we estimated inhomogeneity 
length scale and magnetic field strength of a filament as a whole using period of oscillation and damping 
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time simultaneously. Ruderman & Terradas ( j20~l~3| ) reported that the amplitude decay in standing kink 
oscillations for sufficiently small time is described by a Gaussian function and becomes exponential for 
a later time. The transition from Gaussian to exponential depends on the inhomogeneity length scale. 
Similar studies have been carried out by Pascoe et al. (20121 and Hood et al. ( 2013[ ) for propagating 
waves. Since the knowledge of inhomogeneity length scale is not known a priori thus we use Gaussian 
profile to investigate how it affects the seismology parameters. We find that Alfven speed and magnetic 
field strength are almost same as calculated using exponential damping. However, inhomogeneity length 
scale is different because the Gaussian function falls slowly than exponential function and therefore 
damping time is large. Therefore, we conjecture that based on our MHD seismology that inhomogeneous 
layer was already present when the global kink oscillations are excited. This further causes the quick 
damping of the oscillations via the resonant absorption through the well evolved layer of inhomogeneity 
across the prominence tube. 
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